Cavity cooling of translational and ro- vibrational motion of molecules: 
ab initio-based simulations for OH and NO 
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Abstract We present detailed calculations at the ba- 
sis of our recent proposal for simultaneous cooling the 
rotational, vibrational and external molecular degrees of 
freedom [I]. In this method, the molecular rovibronic 
states are coupled by an intense laser and an optical 
cavity via coherent Raman processes enhanced by the 
strong coupling with the cavity modes. For a prototype 
system, OH, we showed that the translational motion 
is cooled to few /xK and the molecule is brought to the 
internal ground state in about a second. Here, we inves- 
tigate numerically the dependence of the cooling scheme 
on the molecular polarizability, selecting NO as a sec- 
ond example. Furthermore, we demonstrate the general 
applicability of the proposed cooling scheme to initially 
vibrationally and rotationally hot molecular systems. 

PACS 33.80.Ps, 32.80.Lg, 42.50.Pq 



1 Introduction 

Control at the quantum level over both the internal and 
external degrees of freedom of gas-phase molecules has 
been pursued by several groups worldwide [2J. This in- 
tense activity has various ultimate goals, such as the re- 
alization of Bose-Einstein condensation of complex sys- 
tems, the emergence of new ultracold chemistry [3], and 
the applications for high-precision measurements [U__j 
and quantum information processing [6l__J__]. The effi- 
cient preparation of molecules into predetermined quan- 
tum states is the prerequisite for control. Various ap- 
proaches are currently applied. Established methods for 
generating ultracold alkali dimers use photoassociation 
and, or in combination with, Feshbach resonances [5]. 
Another approach is based on buffer gas cooling, in which 
the molecules are thermalized with a cold buffer gas [___] . 
Deceleration of molecules from supersonic nozzles (111 
[12], filtering from an effusive source [13] and collisional 



techniques [TH[T5] are other methods which have been 
developed to produce samples of cold molecules. 

The application of laser cooling techniques, which 
had tremendous impact on atomic and optical physics [161 
fT7lfT8] . to molecules turned out to be inefficient so far. 
Here, multiple scattering channels, due to the molecu- 
lar ro- vibrational structure, are coupled by spontaneous 
emission, leading to uncontrolled heating of the system 
as an undesired side effect. Even when heating can be 
suppressed, the resulting cooling time is long and may 
only be feasible for molecules which are confined in ex- 
ternal traps for very long times |T9][20] . Proposals aimed 
at overcoming this problem have been suggested: Sophis- 
ticated laser-cooling schemes, based on optical pumping 
the ro-vibrational states [2TII22] and optimally modu- 
lated femtosecond excitation pulses [23l[24] , achieve effi- 
ciencies which are still severely limited by spontaneous 
decay and its inherent long time duration. 

An alternative way to optically cooling external de- 
grees of freedom of molecules was proposed and studied 
in [251 26,27.28], and it makes use of the enhancement 
of stimulated photon emission into the cavity mode over 
the spontaneous decay. This mechanism was successfully 
applied for cooling the motion of atoms 29 , 30, 31j . In a 
recent work pQ, we extended this idea and proposed to 
optically cool the external as well as the internal molec- 
ular degrees of freedom by using the coupling to a cav- 
ity. In the proposal the molecules are cooled by a Ra- 
man process in which photons from an intense laser are 
scattered into the modes of a resonator by molecular 
dipolc transition, thereby cooling the motion to the cav- 
ity linewidth and the ro-vibrational degrees of freedom to 
the ground state. The dynamics is numerically simulated 
for the OH radical as a prototype system and for a real- 
istic set of experimental parameters, while the molecular 
properties are calculated with state-of-the-art ab initio 
calculations. Starting with an internal temperature of 
300 K we could show that the translational motion is 
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cooled to a few /xK and the internal state is cooled to its 
ro- vibrational ground state pQ. 

In this article we provide the details of the theory 
and of the ab initio results at the basis of [T] . In order 
to demonstrate explicitly the general applicability of the 
proposed scheme, we present further numerical simula- 
tions where we apply it for the simultaneous cooling of 
vibrational and rotational degrees of freedom for OH in 
a set up with an artificial high temperature, such that 
several vibrational states are initially occupied. More- 
over, we simulate cooling of a similar diatomic system, 
NO, whose 7r bonding allows for a larger electronic po- 
larizability. 

This paper is structured as follows. First the central 
idea is introduced in section [21 followed by a detailed ac- 
count of the theoretical model including the basic equa- 
tions in section [3l In section [4] the two selected mole- 
cules OH and NO are introduced and their properties 
are evaluated with ab initio methods. In section [5] the 
experimental boundary conditions are given which are 
the final ingredients leading to the results presented in 
section [6l The paper closes with a summary, conclusions 
and outlook. 

2 Optically cooling molecules to the ground 
state 

The cooling method we propose is based on the enhance- 
ment of the (anti-Stokes) Raman transitions, where pho- 
tons of an intense laser field are scattered by a molecu- 
lar dipole transition into the modes of a cavity thereby 
cooling the external and the internal molecular motion, 
as sketched in Fig. [TJ Spontaneous emission is largely 
suppressed, since the laser light is far detuned from all 
direct electronic transitions. The enhancement of co- 
herent scattering into the cavity modes is obtained by 



driving the Raman anti-Stokes transitions on resonance. 
For this purpose, the laser frequency is varied sequen- 
tially, in order to shift the comb of cavity resonances 
across the molecular anti-Stokes transitions, as sketched 
in Fig. [2j The external motion is initially cooled to the 
cavity lincwidth by setting the laser frequency on the 
red-side of the Raylcigh line [32J in Fig. [3] In our numer- 
ical simulations we checked that this is a relatively fast 
stage, in which the internal motion is not affected by the 
scattering processes. Then, the ro- vibrational degrees of 
freedom are cooled as the laser frequency is sequentially 
tuned to address all relevant anti-Stokes lines. The ro- 
vibrational levels are hence sequentially emptied, till the 
molecules finally end up in their ground state. 

Rayleigh Peak 
cavity modes 1 ^ ^ 
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Stokes lines t molecular anti-Stokes lines 
laser light 

Fig. 2 Combined with the cavity modes, the laser drives 
a sequence of coherent Raman anti-Stokes transitions which 
on average reduces the ro- vibrational energy of the molecules. 
The spectrum can be shifted relative to the fixed cavity lines 
by tuning the laser frequency. The frequency comb of the 
cavity modes is much denser than in the figure. 



3 Model and basic equations 



a) b) 




Fig. 1 The idea in a nutshell: (a) molecules which are as- 
sumed to be confined in a high-finesse cavity interact with 
strong CW laser light. The laser intensity might be enhanced 
by a build-up cavity, (b) The laser light is far detuned from 
all direct electronic transitions from the ground state X to 
electronically excited states A,B,... in the molecules. The ro- 
vibrational energy of the molecules is reduced by a sequence 
of Raman transition, where laser photons are scattered into 
the cavity modes. 



We now summarize the theoretical model, describing the 
quantum dynamics of the molecular internal degrees of 
freedom and of the modes of the cavity field. We consider 
a gas of molecules of mass M, prepared in the electronic 
ground state X, and with dipole transitions X — > E. 
Here, E = {A, ...} is a set of electronically excited states, 
which may contribute significantly to the total polariz- 
ability. We denote by the ro-vibrational states of 

the electronic ground state, such that the state \X,£) is 
at frequency , and by the ro-vibrational states 
of the electronic excited states, with \E, £") at frequency 
. The elements of the dipole moment d are 

v^, = (E,e\d\x,o. a) 

Denoting by p, x the conjugate momentum and position 
of the center of mass, the Hamiltonian describing the 
dynamics of the uncoupled molecular degrees of freedom 
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reads 



2M 



(2) 



E,£" 

The transitions between molecular levels are driven by 
a far-off resonant laser and interact with an optical res- 
onator as illustrated in Fig. [Th.) , undergoing a coherent 
dynamics which is described by the total Hamiltonian 
H = H mol + H C + V C + V L . Here, 



a\a c 



(3) 



is the Hamiltonian for the cavity modes, where a c (a[) 
denotes the annihilation (creation) of a photon with en- 
ergy hf2 c , with wave vector k c , polarization eo and the 
(vacuum) electric field amplitude £ c . The coupling be- 
tween cavity modes and molecular dipole transitions is 
described by the term 



E h 9 



c,€'-£" 



(z)a c 



(4) 



where <? c .£'— £" (z) = £ c (z)(eo ■ T>^i^^»)/h is the vacuum 
Rabi coupling, whose dependence on the position along 
the cavity z axis is determined by the corresponding spa- 
tial mode function, and 



U" 



(5) 



are the corresponding dipole lowering and rising opera- 
tors. The molecules are pumped by a laser, which cou- 
ples to the dipole and is described by a classical field 
with amplitude £l, polarization cl, frequency lul, and 
wave vector Icl- The laser is assumed to propagate per- 
pendicular to the cavity axis, and the Hamiltonian term 
describing the coupling with the dipole takes the form 

v L = ^L,^e- i( ^ t - kL - x) (^+^ r )+h.c., 

(6) 

where Ql^^i = £l(V^^^> -£l)/^ gives the correspond- 
ing coupling strength. In the interaction Hamiltonian Vl, 
we keep the counterrotating terms, as we will consider 
parameter regimes where the rotating wave approxima- 
tion is not valid. 

Finally, cavity modes and dipole transitions couple 
to the modes of the electromagnetic field external to the 
resonator. We denote by rjt«_^> the decay rate along 
the dipolar transition \E,£") — > \X,£'), such that the 
linewidth of the state \E, £"} is r?, = The 
linewidth of the resonator modes is 2k, and we assume 
the regime in which the Free Spectral Range (FSR) Q c +i- 
Q c > 2k. 

We assume that the laser and the cavity modes are 
far-off resonance from the molecular transitions, as de- 
picted in Fig.QJ)), such that the most relevant processes 



arc Raman transitions, i.e., incoherent Raman scatter- 
ing in which a laser photon is absorbed and then emitted 
spontaneously, and coherent Raman scattering in which 
a laser photon is scattered into one cavity mode. The 
rate of spontaneous Raman scattering reads 

r ^£'(P) = E A"-? (^»(P) + 7^- r (p)) 

(7) 

and it depends on the molecular momentum p through 
the Doppler effect as 



7b*-?' (P) 



Q 2 



(^ £ ,,+k L -p/M)2 + /|„/4 
O 2 



(8) 



(Af £ „ - 2uj l - k L • p/M) 2 + r 2 /4 (9) 



where Af^i 



,(E) 



lol denotes the detuning 



between laser and internal transition. The rate ^jL^ (p) 
gives the scattering rate of a photon from the laser into 
the cavity mode and its subsequent loss from the cavity. 

To evaluate the rate of photon scattering into the 
cavity modes, we consider a standing wave cavity, and fo- 
cus onto the regime in which the cavity photon is not re- 
absorbed by the molecules but lost via cavity decay [26] • 
Moreover, we assume that the molecular kinetic energy 
exceeds the height of the cavity potential. In this regime, 
we decompose the coupling constant <7 Ci £/_>£» (z) into the 
Fourier components g c at cavity- mode wave vector 

±|k c |, and write the rate of photon scattering into the 
cavity as -T^L^, (p) = r£^,(p) + r£^, (p), where the 
sign ± gives the direction of emission along the cavity 
axis and 



T4( P ) 



2. £ 

c,E,£" 



(10) 



7 E ,£-£»(P) 



(<5w+±k c -p/M)' 



± k c • p/M) 2 + K 2 



where 5u> 



XX) 



XX) 



± lol — ft c is the frequency 



difference between initial and final (internal and cavity) 
states. We note that the reabsorption and spontaneous 
emission of the cavity photon can be neglected when 



lsc,e-e"^,«-«"/^fVl 



(ii) 



The basic condition for cavity cooling is that the rate 
of photon scattering into the cavity exceeds the corre- 
sponding spontaneous Raman scattering rate, 

r^ e (p)>>r^,( P ). (12) 

This is the regime on which we will focus. 
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3.1 Rate equations 

We can now write the equations describing the damped 
dynamics of the molecules. For this purpose, we denote 
by Wf (p) the occupation of the molecular state \X, £) at 
momentum p. 

An equation for the kinetic energy can be simply de- 
rived in the scmiclassical limit. Let us first define by 
Ap 2 /2M the average change in kinetic energy in an in- 
finitesimal interval of time At. This is given by Ap 2 /2M — 

Ap\j2M, where 



= At 



dp £ 



(13) 



/(p)W 4 (p) 



and we have introduced, for convenience, the total rate 
rp^ip) = rg'^(p) + rg,^(p)/2. In the regime where 
the laser is tuned such that only Rayleigh scattering into 
the cavity mode is enhanced, the relevant processes will 
not change the internal state of the molecule. In this 
regime, when the rate of spontaneous Raman scattering 
is much smaller than the coherent cavity scattering one 
obtains the typical dynamics of Dopplcr cooling, see [33j 
for a detailed treatment. Here, the cooling rate is given 
by {H 2 k 2 /2M)r'(0) 1 where f (0) « ^ 1^(0) and the 
asymptotic value of the kinetic energy is 



= h {uj l - n c ) 2 + k 2 
kin 4 \u L -n c \ ' 



(14) 



which is minimum for u>l — fl c = —k and takes the value 

EST* = Hk/2 m- 

The generic rate equation, describing the dynamics 
of the population of the molecular state V$ = J dpW^ (p) 
reads 

^ = - /*£ fere?) + r l^(p)) WC(P) 
+ [ d P J2 (r^{p) + ^(p)) w e (p), 

(15) 

which requires the knowledge of all rates and popula- 
tions at each instant of time. Assuming that the mole- 
cules have been previously cavity cooled to the cavity 
lincwidth, and that the momentum distribution is not 
significantly modified by the inelastic processes that change 
the internal state of the molecule, then the rate of inelas- 
tic processes is 



r;4,(o)% 



Hence, in this regime, the incoherent dynamics can be 
described by the rate equation 



r 7 



r 7 



(17) 



(16) 



whose stationary solution is found using the detailed bal- 
ance principle. 

The simulation of the cooling procedure uses rate 
equations for the molecular levels with the rates ([7]) 
and (flU)) as evaluated in the previous section. The se- 
quence for the laser driven transitions is simulated by 
changing the detuning at the different steps. The cool- 
ing of the external degrees of freedom is the preliminary 
step. We have evaluated the cooling rate for the external 
motion with the laser frequency set to lol — f2 c — k/2, 
and checked that the internal molecular level distribu- 
tion is practically unaffected during the cooling time. In 
the step for cooling the internal degrees of freedom, we 
then assumed that the motion was cooled at the Dopplcr 
limit, and set p = in the rates ||7J) and (fTU|) . In the next 
section we report the details of the numerical treatment, 
for determining the parameters which give us the rates. 



4 The molecular systems: OH and NO 

We simulate the cooling scheme for OH and NO radicals. 
Cold ensembles of these two molecules are experimen- 
tally available [T4 1 I34 1 I35 1 I36] . We first apply the cooling 
scheme for an internal temperature of 300 K. Since at 
this temperature OH molecules are in the vibrational 
ground state, we also simulate the process for OH in a 
situation in which initially the internal temperature is 
artificially high, in order to investigate cooling of vibra- 
tional motion explicitly. 

OH and NO can be modeled accurately by ab initio 
quantum chemistry methods. They have a similar elec- 
tronic valence structure which allows to compare their 
electronic properties. In particular, the unpaired electron 
leads to a significant spin orbit splitting of 139 cm -1 and 
123 cm" 1 for the X 2 II ground states of OH and NO, 
respectively, yl-doubling occurs due to the coupling be- 
tween rotation and electronic angular momentum. The 
doubling lies in the MHz regime and is different for 
the X 2 7T 3 /2 and the X 2 II 1 i2 components. We only con- 
sider one yl-state, since states of different parities are not 
coupled by Raman transitions. The hypcrfinc splittings 
in OH, due to the nuclear spin of the hydrogen atom, 
are neglected since angular momentum conservation for 
the rotational Raman transitions inhibits transitions be- 
tween these sublevels. Hence, we assume that it is possi- 
ble to prepare the molecules in the lower lying X 2 n 3 / 2 
state for OH and X 2 n 1 / 2 for NO. With these assump- 
tions we choose J as the rotational quantum number. 
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Because of the Raman selection rule A J = 0, ±2, matrix 
elements between rotational states with even and odd J 
values vanish. This leads to two separate ladders for the 
scattering process with ground states \X, v = 0, J = 0) 
and \X, v = 0, J = 1), respectively. 

The energy of the rotational levels reads 



Ej = B(v)J{J + 1) - D r J 2 {J + l) 2 



with 



B(v) = B e 



B e . x2 {v+-f 



(18) 



(19) 



The constants, including the anharmonicitics, are taken 
from Ref. [37]. In particular, for OH they read(i? e , B B}X i, 
B e , x2 , toe) = (18.871, -0.714, 0.0035, 3735.21) cm" 1 while 
for NO they are (B e , B e>xl , B BtX2 , u e ) = (1.7042, -0.01728, 
0.000037, 1904.2) cm' 1 . The anharmonicity in the vibra- 
tional ladder is included in the ab initio potential energy 
surfaces (PES). 



in which the integration is over the nuclear coordinates R 
and where indicates the pure vibrational states. The 
transition strength between the rotational levels is evalu- 
ated using the Placzek- Teller coefficients [3"51l3"9~] for gen- 
eral diatomic molecules (shown in Tab. [I]). The Placzek- 
Teller coefficient Sj—>ji for the corresponding transition 
is multiplied with the square of the absolute value of 
the polarizability, giving the complete transition matrix 
element 

\a v ,j^v',J'\ = \a v ^ v '\ 2 ■ Sj-,j>. (24) 

The obtained value for a can be inserted in Eq. [7] and 
Eq. [TO] to calculate the transition rates in the molecular 
system. 



Table 1 Placzek- Teller coefficients S adapted to diatomic 
molecules. With the different cases for J one can calculate 
Stokes and anti-Stokes transitions as well as Rayleigh scat- 
tering. 



4-1 Molecular polarizabilities 

We now determine the molecular rotational and vibra- 
tional structure in order to evaluate the scattering rates 
for simulating the cooling dynamics. A relevant quantity 
is the polarizability a, which is defined for molecular sys- 
tems as 



E 

E,t" 



V. 



V. 



•Vf 



h{A E 



(20) 



AJ 



J + 2 



J — 2 



S 



[J(J+1) 



J(J + l)(2J-l)(2J + 3) 

3(J + l) 2 (J + 2) 2 
2(J + l)(J + 2)(2J + l)(2J + 3) 

3(J- if J 2 
2(J- 1)J(2J + 1)(2J- 1) 



and is expressed as a function of the laser frequency ujl 
and the internal nuclear coordinates. Note that also higcr 
lying excited states can have significant contributions. 
The polarizability is described as a tensor which can 
be expressed in terms of irreducible components. The 
isotropic part , defined as 



a (0) = brr(a), 



(21) 



is responsible for transitions which do not change the 
angular momentum of the molecule. The tracelcss com- 

(2) 

ponent a zz along the main axis of the molecule, can be 
evaluated from 

= afi = 



\ ( a ij 



(22) 



This component is needed to describe ro-vibrational tran- 
sitions associated with a change in J. In all subsequent 
formulas a represents the appropriate component. This 
general formulation allows us to treat all possible tran- 
sitions between the internal states. 

The transitions between the vibrational states can be 
written as the matrix clement 



Rayleigh 

anti-Stokes 




(Z v \a{R)\&), 



(23) 



Fig. 3 Illustration of the Placzek- Teller coefficients 5 for 
the seven relevant Raman lines. The transitions indexed with 
n, which are allowed according to the selection rule A J = 
±2 are shown. For example, —AJn — 2 refers to J2^o and 
— AJn — 4 refers to J^^i. The coefficients for the Rayleigh 
scattering AJ — are not included in the picture. 



4-2 ab initio calculations 

We carried out different types of ab initio calculations 
for all the relevant properties not directly available from 
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6 8 10 12 14 16 
v [GHz] 

Fig. 4 Reduced spectrum of the seven relevant Raman lines 
for OH. The Raman transitions are projected into one FSR. 
The solid lines indicate the driven anti-Stokes transitions. 
The dashed lines represent the corresponding Stokes lines in 
the Spectrum. The cavity resonances are the lines on left and 
on the right side highlighted by the arrows. 

literature. The PES were obtained with Molpro [40] in 
a complete active space self consistent field (CASSCF) 
calculation and were further improved with a multi ref- 
erence configuration interaction (MRCI) step. A valence 
triple zcta single atom basis set (cc-pVTZ) [5T] was used. 
For the OH radical we chose an active space of 11 or- 
bitals which includes all electrons except the Is core or- 
bital from the oxygen. The NO was treated in a similar 
way with 12 active orbitals and only leaving out the Is 
core orbitals from the oxygen and the nitrogen atom. 

For the calculation of the polarizability linear re- 
sponse theory on the same CASSCF [421I43] level were 
performed. This method is implemented in the program 
package DALTON [44] . The quality and convergence of 
the resulting polarizabilities were checked through sev- 
eral tests. The finite field method is used for comparison 
with the values in the static limit where ujl = 0. The 
static polarizability of a particle can expanded in a Tay- 
lor series with respect to the perturbation of an electric 
field, such that a is found from 



, FF 



d 2 E(£) 
d£id£j 



(25) 



The differentiation is done numerically with small field 
strengths. With this method several basis set checks (see 
Tab. © were performed for OH on the CASSCF level 
mentioned before. It turned out (as already mentioned 
in earlier publications [46]) that the results are very sen- 
sitive to the choice of the basis set. The addition of dif- 
fuse functions gives the better results. Higher angular 
momentum f-polarization functions have only a little ef- 
fect on the polarizability. It turned out that correlation 
consistent basis sets (cc-pVXZ) [4T1I47] are more ade- 
quate for this task than the widely used 6-31G basis sets 
[4"5] , The aug-cc-pVTZ basis set was chosen for the po- 
larizability calculations as a good compromise between 
size and quality of the result. This basis set was used 
to compare the result from the finite field calculations 
against linear response calculation in the static limit at 



Table 2 Single point polarizability calculations for OH with 
various basis sets. All calculations were carried out with on 
the CASSCF level of theory with the finite field method 
(£ — 0.001 au) at the equilibrium geometry. Listed are the 
needed tensor components a' ' and a^J . Additionally the 
corresponding dipole moment is printed (The experimental 
value is /i z = 0.657 au= 1.67Debye [H]). All values are given 
in atomic units. 



Basis set 




(2) 




6-31G(d,p) 


3.809 


2.172 


0.711 


6-31G++(d,p) 


4.936 


1.652 


0.741 


aug-cc-pVDZ 


6.826 


1.438 


0.648 


cc-pVTZ(d,p) 


5.194 


1.836 


0.667 


cc-pVTZ 


5.105 


1.928 


0.668 


aug-cc-pVTZ(d,p) 


7.230 


1.188 


0.648 


aug-cc-pVTZ 


7.240 


1.200 


0.645 


aug-cc-pVQZ 


7.407 


1.069 


0.646 


aug-cc-pVQZ(d,p) 


7.403 


1.060 


0.647 



Table 3 Comparison between the OH polarizabilities from 
the finite field and the linear response calculations at the 
CASSCF/aug-cc-pVTZ(d,p) level of theory. Additionally the 
dynamic polarizability at 532 nm is shown. All values are 
given in atomic units. 



Method 




a {2) 


Finite field 




7.23 


1.19 


Linear response, 


static limit 


7.21 


1.17 


Linear response, 


A = 532 nm 


7.39 


1.15 



16 



12 - 



Z5 
CO 













a (0) — 

(2) 

cc 1 — - 


0.6 0.8 


1 1.2 1.4 1.6 1.8 2 2.2 



r(OH) [A] 

Fig. 5 The polarizability function in dependence of the bond 
length for OH at a laser wavelength of 532 nm. The isotropic 
part a' ' is used for the calculation of the Rayleigh scattering 

(2) 

transitions. The traceless part a z J of the tensor determines 
the rotational transitions. The equilibrium geometry is indi- 
cated by tq. The wiggles between 1.6 A and 2.0 A are artifacts 
from resonances where the detuning Af^i is near zero. These 
singularities were smoothed by an interpolation. 



zero frequency, where both methods should give equal 
results (Tab. [3]). The obtained values differ less than 1 % 
and hence we can assume to have reasonable values. 
The calculated polarizability function in dependence of 
the OH bond length is shown in Fig. [5] The function 
has been calculated for a laser wavelength of 532 nm. 
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Both curves for the isotropic and the traceless part show 
a typical increase with greater bond lengths reaching a 
maximum between 1.6 A and 2.0 A. In this area the X 2 U 
and A 2 E + state come closer (see Fig.Q]) and their energy 
difference matches the laser frequency. Here, the Raman 
approximation |49j becomes invalid and the evaluated 
data points show a singularity. This effect is smoothed 
by an interpolation routine. The remaining wiggles have 
no effect on the results as they lie outside the nuclear 
range relevant in the considered Raman process. 

The vibrational eigenvalues and eigenfunctions of the 
molecules were evaluated with a relaxation method us- 
ing propagation in imaginary time with an additional 
diagonalization step [50]. This method is used to deter- 
mine the vibrational wavefunctions in the ground state 
and first electronic excited state, necessary to obtain the 
relevant Franck-Condon factor of the system. 

5 Experimental parameters 

We now discuss the experimental parameters, which are 
required in order to implement efficiently the cooling 
scheme. We assume that the OH or NO molecules are 
prepared in the lower- lying component of the X 2 II elec- 
tronic ground state with a motional temperature around 
or slightly below IK. This could be realized with, e.g., 
helium-buffer-gas cooling [10] , electrostatic filtering [13] , 
or decelerator techniques [TT|I12|. However, as the pur- 
pose of this paper is to demonstrate in particular the 
cooling of the internal motion of the molecule, the filter- 
ing method [13] matches most closely our assumed initial 
conditions. At T « 300 K only the first nine rotational 
levels of OH are occupied with a population > 10~ 3 . 
The first vibronic excitation for OH has an energy of 
3727.95cm -1 [37], and thus only the ground state v = 
is occupied. For NO there are 25 rotational levels and be- 
cause the vibronic excitation lies at 1904.2 cm -1 again 
only the ground state is occupied. 

We choose a laser wavelength of 532 nm. Since this 
can be generated by frequency doubling of 1064 nm light, 
commercial systems are available with ample power, with 
high quality of the beam profile and which run at single 
frequency. The frequency of this light is far below that of 
the OH A-X ro- vibronic band. We assume to have single- 
frequency light of 10W enhanced by a factor of 100 by 
a build-up cavity in a TEMoo mode, corresponding to 
a Rabi coupling £2i,fi— >o — 27r x 69 GHz and frequency 
uj l = J A) - J X) - A with A ^ 2n x 407 THz. Note 
that these intracavity intensities arc routinely achieved 
in some commercial devices. 

The laser frequency should be varied during cooling, 
in order to drive (quasi-) resonantly the cooling transi- 
tions. In combination with the broad spectrum of cav- 
ity modes, the laser only needs to be varied over one 
FSR to address all anti-Stokes lines. Fig. 2] displays the 
Stokes and anti-Stokes Raman lines as a function of the 



frequency modulus the FSR. Most commercial systems 
delivering high-power 532 nm light are not tunable, but 
this does not seem a fundamental problem. The build- 
up cavity will require tracking while the pump laser is 
scanned, which is also possible with today's technology. 

We choose a high-fincsse optical Fabry-Perot-type 
cavity with a length of L = 1 cm, which seems large 
enough to accommodate a trapping device as e.g., demon- 
strated in refs. [T9 l l20 p 5i p 52j. This length leads to a free- 
spectral range (FSR) of 15 x 27rGHz. For simplicity we 
assume that the cavity only supports zeroth-order trans- 
verse modes. The cavity half-linewidth is set to K = 
75 x 27rkHz, and the coupling <? C) o— >o = 27r x 116 kHz. 
This is achieved with a mode volume of 3.2 x 10 _13 m 3 , 
assuming a mode waist of wq = 6 ^m, and a cavity finesse 
F = 10 5 , i.e., a mirror reflectivity of 0.999969. In practice 
it might be necessary to choose a near-degenerate setup, 
e.g. a near-confocal cavity, in order to avoid higher-order 
cavity modes from appearing everywhere in the reduced 
spectrum, possibly coinciding with Stokes lines. Choos- 
ing such a degenerate setup actually increases the pos- 
sibilities of scattering light into the cavity, increasing 
the cooling rate. There is one caveat: if the higher-order 
transverse cavity modes are not exactly degenerate with 
the zeroth order, this will lead to a cavity linewidth 
which will appear broadened and all lines in the re- 
duced spectrum should be convoluted with this cavity 
linewidth. However, we do not expect this to be an is- 
sue, not even in a near-degenerate cavity, as our scheme 
is robust against a small number of coincidences between 
Stokes and anti-Stokes lines. 



6 Cooling simulation results 

6. 1 Rotational cooling of OH and NO radicals 

We now simulate the cooling of OH and NO with the 
experimental parameters as described above. We start 
with the OH radical as the molecule of main interest. 
Fig. [S] and Fig. [5] show the results of the simulation. The 
FSR has been fine tuned so that the transitions J3— >i 
and J2^q could be pumped simultaneously. In general, 
the spectral lines of a molecular system are not equally 
spaced and in general only one transition at a time can 
be addressed. However, two transitions can be selected 
simultaneously by tuning the FSR so that two cavity 
lines coincide with two transitions. 

In Fig. [6] two different sequences are shown. The 
dashed curve represents a simple "top-down" variant 
which cools down the molecule to the ground states be- 
ginning from the highest occupied level. The solid curve 
is a manually optimized sequence in order to obtain a 
faster decrease of the expectation value (J) PjJ. 

Snapshots of the population time evolution are shown 
in Fig. [7] and Fig. [5] In the top-down sequence, Fig. 
[3 we start driving the transitions corresponding to the 
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Optimized sequence 
Top-down sequence 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
time [s] 

Fig. 6 Rotational cooling of OH starting from 300 K. The 
expectation value (J) = Y2jVjJ of the rotational quantum 
number is plotted against time. The solid line shows a cool- 
ing process with a manually optimized laser sequence. The 
dashed line represents a "top-down" process where all higher 
levels are emptied sequentially in multiple cycles. 
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Fig. 8 Populations at selected times of the cooling pro- 
cess for the manually optimized sequence. Each transition 
is driven for 60 ms. 



highest occupied states, Js—>e and J7—>5, continuing with 
the subsequent anti-Stokes transitions. Every step lasts 
about 60 ms, whereby this optimal value was found em- 
pirically. After the first cycle 88 % of the population is 
in the states J = 0,1. To achieve a lower final tem- 
perature the procedure is repeated six times. The rate 
for the decrease in (J) is 2.3 Hz. In the optimized se- 
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Fig. 7 Populations at each step of the cooling process for 
the "top-down" sequence in the first cycle. Each transition is 
driven for 60 ms. 

quence, first the J = 2 and J = 3 levels are addressed. 
In Fig. [5] this can be seen in the snapshot at 120 ms. 
Here a large population fraction is transferred to the 
two ground states. In the next step (after 240 ms) the 
transitions Js^3 and J±—,2 are driven sequentially for 
60 ms each. Repeating the first step, more than 90% of 
the population is transferred to the ground states. By 
addressing the levels J > 5 and repeating the whole 
procedure, 98.8 % of the population is transferred to the 
two ground states J = and J = 1 after 1680 ms. The 
achieved rate for the decrease in (J) is 3.6 Hz. 

We also checked the effect of spontaneous Raman- 
processes on the vibrational state distribution when the 
rotational motion is cooled, and found only little vibra- 
tional heating. After a complete sequence of 1.8 s less 



than 1 % is in the state v = 1. Without cavity cool- 
ing, the rotations are heated, such that, the expectation 
value (J) rises from 2.44 to 2.56 on time scales of the 
order of 1.8 s. 

To analyze the dependence of the cooling process on 
the related molecular properties, we now simulate cool- 
ing of a NO molecule, whose level structure in the elec- 
tronic ground state is similar to the OH radical while the 
rotational and vibrational constants differ. In contrast to 
the more elaborate calculations of OH, the polarizability 
of NO was calculated in the quasi-static limit. This is a 
reasonable approximation as long as the ground state is 
coupled far-off resonance from the excited state. Due to 
the 7r-bonding in the NO radical the value for the polar- 
izability is about 4 times greater than in the OH radi- 
cal. Hence, compared to OH, NO offers a larger polar- 
izability, which increases the cooling rate. On the other 
hand, due to the smaller rotational constant 25 rota- 
tional states are occupied at 300 K. Thus, NO is charac- 
terized by a broader distribution over the rotational lev- 
els than OH, which increases the cooling time. The result 
for the cooling procedure is shown in Fig. [9] and Fig. [10] 
for the manually-optimized sequence. It turned out that 
it is a good choice to start the sequence by addressing the 
levels near the maximum of the initial Boltzmann distri- 
bution (Fig. ITU)) , in this case at levels J = 12, 13. Every 
transition is driven for 5 ms. After 60 ms the maximum of 
the distribution function is shifted to the ground-state 
levels J = 0,1. The next step is started in the states 
J = 20, 19 and after 155 ms more than 50% are in the 
two defined ground states. After 965 ms 99.2% of the 
population is in the lowest two states. The optimal rate 
for the decrease in (J) is 5.3 Hz. In order to make a 
comparison with OH, we repeat the simulation for NO 
by taking the same initial internal distribution of only 
nine levels (which would correspond, in an experiment, 
to a sample of NO molecules at about 25 K). Due to the 
larger polarizability of NO, the cooling is about 12 times 
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0.2 0.4 0.6 0.8 1 



time [s] 

Fig. 9 Rotational cooling of NO with an initial tempera- 
ture of 300 K. The sequence for driving the transitions was 
manually optimized. Due to the larger polarizability the NO 
radical can be cooled faster than OH although more states 
are occupied. 




0.2 0.4 0.6 0.8 1 1.2 1.4 
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Fig. 11 Cooling of the vibrational motion of OH for an ini- 
tial thermal distribution where the first nine levels are sig- 
nificantly occupied. 
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Fig. 10 Populations during the cooling process for NO with 
an optimized sequence. The initial distribution corresponds 
to a temperature of 300 K. Various stages of the cooling pro- 
cess are shown. 



faster than in the case of OH, and the optimal rate for 
the decrease in (J) is 43 Hz. 



6.2 Cooling of the vibrational and rotational motion 

To test the efficiency of our scheme for cooling vibra- 
tional excitations, we consider again OH molecules and 
construct a scenario where the vibrational temperature 
is artifically high, taking an initial distribution in which 
the first nine vibrational levels are appreciably occupied. 
First we just cool the vibrational degrees of freedom, ne- 
glecting the Placzek- Teller coefficients. In this case, not 
only the absolute value of the polarizability, but also its 
finite slope with respect to the internuclear coordinate 
(see Fig(5|) is important. 

A roughly optimized frequency sequence for this pro- 
cedure leads to the fast cooling process shown in Fig.fTTl 
demonstrating that the concept works. However, it over- 
estimates the cooling rate because the neglected Placzek- 
Teller coefficients for anti-Stokes transitions are smaller 
than 0.5. 



In the next step the rotational degrees of freedom 
are combined with two vibrational levels each contain- 
ing nine rotational substates. All initial population is in 
v = 1, with a 300 K rotational temperature. The task of 
achieving a good result for both {v ) and (J) is now more 
complicated but still feasible. The roughly optimized se- 
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Fig. 12 The expectation value of n as function of time, 
where n is a substitute for the quantum numbers v and J, 
respectively. All population is initially in vibrational state 
v = 1 and the rotational levels follow a thermal distribution 
at 300 K. 



quence of Fig |12l is an example where vibrational and 
rotational degrees of freedom have been cooled simul- 
taneously. For these results, the selection rules are ex- 
ploited such as to reduce (v) and (J) simultaneously as 
much as possible. However, as can be seen in the figure, 
this is not always possible. In the beginning the tran- 
sitions vi-^qJ^i and vi_»oJ2->o are driven for 100 ms 
each. In the next step the transition vi^q Ji^i is driven. 
Here (v) still becomes smaller but (J) stays constant. In 
the following step after, 300 ms, a compromise was made. 
The state v = 1, J = has to be emptied but the logical 
target state v = 0, J = is not accessible because the 
Placzek- Teller coefficient for the Jq^o transition is zero. 
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Therefore the next accessible level, v = 0, J = 2, is used, 
which results in an increase in (J). Note that J changes 
by 2 and v by 1 and thus the variations in (J) are larger 
than those in (v). The further progression for the laser 
induced transitions implies the alternation of two prin- 
ciples: transferring populations to the lower vibrational 
states like for the first cooling steps and cooling the ro- 
tational levels without vibrational heating. It turned out 
that for the further cooling steps it is more effective to 
first cool the rotational motion without changing the vi- 
brational motion before transferring the population to 
the vibrational ground state. 

In the end of the procedure (v) and (J) converge to 
their respective minima and 0.5. The achieved rates for 
teh decrease of the expectation values are r v = 0.76 Hz 
and Tj = 0.49 Hz. 

7 Conclusions and outlook 

We have extensively discussed a scheme which allows 
to cool simultaneously the internal and external motion 
of molecules. In particular, we have detailed the rate 
equations at the basis of the simulations and the ab 
initio calculations of the molecular properties. Special 
care was taken to use a realistic set of parameters. For 
these parameters, cooling of the external motion from 
1 K to a few \i K can be accomplished in about a sec- 
ond, thus requiring the support of trapping technologies 
which are stable over these times [TUlH^l^IIl[5T l[52l[55] . 
We performed extensive calculations on the OH radical, 
and studied the cooling efficiency on the molecular prop- 
erties by comparing the cooling efficiency of OH with 
the one of NO. We see that, even if NO has more popu- 
lated rotational levels than OH at the same initial tem- 
perature, NO can be cooled faster because of its higher 
polarizability. We notice that two properties of the po- 
larizability play a different role in the cooling dynamics. 
In particular, the absolute value of the polarizability is 
responsible for the processes which cool the rotational 
excitations, while for vibrational cooling a finite slope of 
the polarizability as a function of the nuclear coordinate 
is essential. 

An interesting question is how the method performs 
for more complex molecules. An advantage of larger mole- 
cules is that states of even and odd J can be coupled by 
Raman processes, contrary to the case of linear mole- 
cules such as OH and NO. On the other hand, with 
increasing complexity, the heating due to spontaneous 
Raman scattering will become more important and must 
be compensated by a faster coupling to the environment 
with the help of a better cavity. Additionally, the cooling 
process can be improved by optimizing the cooling se- 
quence, a task which becomes increasingly difficult with 
more complex molecules and might require advanced op- 
timization strategies like genetic algorithms [M]. Alter- 
natively, our technique could be extended by combining 



vacuum-stimulated emission into the cavity with suit- 
ably designed Raman pulses using optimal control tech- 
niques [23j . 
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